Discrimination of unitary transformations in the Deutsch-Jozsa algorithm. 
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We describe a general framework for regarding oracle-assisted quantum algorithms as tools for 
discriminating between unitary transformations. We apply this to the Deutsch-Jozsa problem and 
derive all possible quantum algorithms which solve the problem with certainty using oracle unitaries 
in a particular form. We also use this to show that any quantum algorithm that solves the Deutsch- 
Jozsa problem starting with a quantum system in a particular class of initial, thermal equilibrium- 
based states of the type encountered in solution state NMR can only succeed with greater probability 
than a classical algorithm when the problem size exceeds n ~ 10°. 
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I. INTRODUCTION 

Quantum algorithms [l[ are typically described in 
terms of the evolution of the state of a quantum sys- 
tem under a prescribed sequence of unitary transforma- 
tions, followed by the extraction of a problem solution 
from the outcome of a measurement performed on the 
quantum system. While the unitary transformations in- 
volved are clearly crucial, they are not the primary ref- 
erence point of the analysis in such a circuit formulation 
of quantum algorithms; this belongs to the state of the 
system. However, in the circuit formulation of any par- 
ticular quantum algorithm, the ingredient which varies 
from one application of the algorithm to another (i.e. the 
input) is typically a unitary transformation. Therefore, 
it seems appropriate to focus on the unitaries in a quan- 
tum algorithm and to regard the algorithm as a tool for 
discrimination of unitary transformations. 

This approach has been suggested before @, Q and 
subject to analysis in selected scenarios 0, @. The pur- 
pose of the latter articles was to extend the Deutsch- 
Jozsa algorithm by investigating the possibility of dis- 
criminating amongst a larger class of unitary transforma- 
tions than that encountered in the original Deutsch-Jozsa 
algorithm. Our purpose is to further promote unitary 
(and state) discrimination as tools for analyzing quan- 
tum algorithms. Specifically we consider the Deutsch- 
Jozsa problem, which is to be solved with the aid of an 
oracle unitary of a specific form. We ask whether the 
notions of unitary discrimination can be used to reach 
the standard quantum algorithm for solving the problem, 
whether they yield alternative algorithms for solving the 
problem and what restrictions they impose on quantum 
algorithms which are implemented on quantum systems 
initially in noisy mixed states. 

The remainder of this paper is organized as follows. 
In Sec. |TT] we describe how oracle-assisted quantum al- 
gorithms can be viewed as unitary discrimination tools. 
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This is applied to the Deutsch-Jozsa problem, using a 
particular oracle unitary, in two distinct ways in the fol- 
lowing sections. In Sec. IIIII wc use unitary and state 
discrimination to arrive at the set of all algorithms 
which solve the Deutsch-Jozsa problem with certainty. 
In Sec. IIVI we assume a restricted set of possible initial 
states, which are mixed, for solving the Deutsch-Jozsa 
problem. We determine a lower bound on the problem 
size, beneath which a classical algorithm will succeed in 
solving the Deutsch-Jozsa problem with greater certainty 
than any quantum algorithm. 



II. DISCRIMINATION OF QUANTUM 
OPERATIONS IN ORACLE ALGORITHMS 

Certain computational problems, such as searching, Si- 
mon's problem Q and the Deutsch-Jozsa problem [J, are 
oracle-assisted, meaning that they are to be solved using 
a binary oracle function / : {0,1}" n- {0,l} m , whose 
form depends on the nature (and the specific instance 
that is invoked) of the computational problem. The task 
is to solve the problem with the fewest oracle invocations; 
the efficiency of the solution is quantified by the number 
of oracle invocations used. Different instances of a given 
problem correspond to different oracle functions. 

The associated quantum algorithms (6 9] require a 
well-determined number, depending on n, of qubits pre- 
pared in a suitable initial state \^o)- These qubits are 
made to evolve collectively in way described by a spe- 
cific sequence of unitary transformations. Ultimately a 
measurement yields an outcome from which the prob- 
lem solution can be extracted with high probability. In 
oracle-based quantum algorithms the oracle is invoked 
via evolution described by a unitary transformation Uf, 
whose form depends on the problem as well as the par- 
ticular oracle function, /, in use, i.e. the particular in- 
stance of the problem. In the simplest cases, to which 
we restrict our consideration, the target of the oracle is a 
single binary variable, i.e. m = 1. Examples include the 
Deutsch-Jozsa 0] and Grover's search algorithm For 
a given type of problem, there can be different possibili- 
ties for the number of qubits required and the structure 
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of the oracle unitary. We consider instances where the 
number of qubits required is n. In terms of computational 
basis states, \x) := \x n ) <£>...&> |xi) with Xi € {0, 1} we 
assume that the oracle unitary operates via 

U f \x) = \x) (1) 

and this is extended linearly to all linear combinations of 
computational basis states. 

Applications of the oracle unitary may be interspersed 
with other unitaries, Vq, . . . , Vm which arc oracle inde- 
pendent, i.e. these remain fixed for all choices of ora- 
cle function /. The resulting algorithm has the struc- 
ture illustrated in Fig. [T] and the algorithm unitary is 
U a \ g = VitfUf . . .UfViUfVo. This is the most general 
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FIG. 1. General structure for oracle-based quantum algo- 
rithm where the oracle output is a single bit. The oracle is 
invoked M times. The initial and final unitaries, Vo, Vm are 
not generally necessary but are enable initialization and mea- 
surement in the computational basis. 

structure for m = 1 oracle-assisted algorithms since two 
successive applications of the oracle unitary provide a 
trivial identity unitary. The final state \^ f) = t/aigl^o) 
depends on the particular oracle function /. 

In such scenarios the input to the algorithm is the or- 
acle function and not the initial state. The algorithm 
output then identifies or classifies the input oracle func- 
tion. For example, in the Dcutsch-Jozsa problem the al- 
gorithm determines whether the oracle function belongs 
to the class of constant or balanced functions (these terms 
are described below). In the algorithm for searching a 
database with one marked item, located at s, the oracle 
is defined as f(x) = if x ^ s and f(s) = 1. Determin- 
ing s is equivalent to identifying which of the possible 
oracle functions was used. Analogously, the associated 
quantum algorithms amount to tools for classifying or 
discriminating between the possible oracle unitaries. 

The framework for unitary discrimination requires a 

set of known unitaries |i7i, U2, ■ ■ ■ j arid associated prob- 
abilities {pi,p2, . . .} . One party chooses one of these uni- 
taries, Uj, with probability pj and another party must 
determine which unitary was selected by applying the 
unitary one or more times, together with other quantum 
operations, to a quantum system. Ultimately a mea- 
surement is performed on the quantum system and the 
choice of unitary is inferred from the measurement out- 
come. Unitary discrimination can be reduced to a quan- 
tum state discrimination problem [jj Hol - [l^ | by apply- 
ing the unitary to a standard initial state, most gen- 



erally described by density operator po, and attempt- 
ing to discriminate between the possible resulting out- 
put states. In this article we consider cases where the 
unitary is applied once only. Thus the possible output 
states are pj = UjPoUj , occurring with probability p 3 -, for 
j = 1,2,.... The general framework for discriminating 
between states [H|,[ll[ requires a POVM with positive op- 
erator elements {tci, 7T2, . . .} that satisfy • Uj = I and 
a rule for associating states with measurement outcomes. 
In the minimum error discrimination scenario, which we 
consider here, we arc required to select one state for each 
outcome and the "undecided" inference of the unambigu- 
ous discrimination (h"H scenario is not permitted. Here 
it is possible to make an incorrect inference and the task 
is to choose measurement that minimize the probability 
with which such an error occurs. For unitary discrimi- 
nation, both the measurement and the initial state must 
be chosen so as to minimize the error probability. 

III. APPLICATION TO THE DEUTSCH-JOZSA 
ALGORITHM ON ARBITRARY INITIAL STATES 

The Deutsch-Jozsa problem consides oracle functions 
/ : {0,1}™ 1 — ^ {0,1} that are required to fall into one of 
two classes: constant, meaning that / returns the same 
value for all possible arguments or balanced, meaning that 
/ returns for exactly half of the arguments and 1 for 
the other half. The task is to determine the function 
class with the minimum number of oracle invocations. A 
classical algorithm proceeds by evaluating / at randomly 
chosen distinct arguments. This will succeed with cer- 
tainty in all cases after N/2 + 1 oracle invocations [7lll5j|. 
where N = 2™ is the number of possible argument val- 
ues. A quantum algorithm exists and @, HH, m its mod- 
ified form [la ], uses an oracle unitary of the form given 
in Eq. ( ED) e xactly once to solve the problem with cer- 
tainty [7|, [l5| , giving an exponential speed-up in terms of 
n. 

We aim to use the notions of unitary discrimination 
to arrive at all quantum algorithms which can determine 
function class with certainty while using n qubits and an 
oracle of the form of Eq. ([1]). Since each function class 
is represented by many unitaries, this requires discrimi- 
nation between two quantum operations. For each given 
class of functions, the quantum operation is 

Po^Pf-= PfUfpoUj (2) 

/ in class 

where the sum is over all possible functions in the given 
class and pf is the probability with which each function 
could be chosen given that the particular class is chosen. 
For constant functions, Eq. (p} implies that Ut = I and 
thus 

PO <-> Pf const = Po- (3) 

For balanced functions, expanding in the computational 
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basis, po = Y.x V toPoxy\x){y\ 1 gives 



N-l 



/3o^Pfbai= E Pf E (-l) /(x)+/(y Vox,|x)(y|. 

/balanced a:,y— 

(4) 

which gives the density matrix elements after the bal- 
anced function operation pfbaixj/ := (x|pfbai|y) as 



Pfbalzy^ E 

/ balanced 



f(x) + f(v) 



PO : 



(5) 



Thus pfbaix* = Po xx for a; = 0, . . . , N - 1. For non- 
diagonal density matrix entries, the summation will be 
complicated by the possibly different probabilities for 
each balanced function. We shall assume that the prob- 
abilities with which each balanced function is selected 
are identical. Thus Pf = 1/B where B is the number 
of balance functions. Enumerating the number of bal- 
anced function amounts to counting the number of ways 
in which N/2 of the N possible arguments which will 
return can be selected. Thus B = (^ 2 ) and 



Pf bal xy — -q PO xy 



E < 

/ balanced 



-i) 



(6) 



The sum £/ balanced (-l) f{x]+f{v) can be evaluated for 
given values of x and y by determining for how many 
balanced functions f(x) = f(y) and for how many f{x) ^ 
f(y). If x ^ y, this is independent of the choices of x and 
y. This can be established by noting that the collection 
of all balanced functions can be listed by assigning to 
N/2 of the N possible argument "slots" and 1 to the 
remaining slots. The argument values, x = 0, 1, 2, . . . N— 
1 in this process merely serve as labels and interchanging 
two of them will not affect the sum. Thus it suffices to 
compute this for x = and x = 1. The four possibilities 
are tabulated in Table [TJ 



f(x = 0) f(y = 1) (-1) 



Number of Instances 



1 

-1 
-1 
1 



(N-2\ 
\ N/2 J 

\N/2-l) 

( N ~ 2 1 
\N/2-l) 

/AT-2N 

V N/2 ) 



TABLE I. Four possibilities for {-l) f ^+l M for x = and 
y = 1. The two leftmost columns provide the possible com- 
binations of values returned by /. The last column lists the 
number of times that each possibility occurs. For example, 
the balanced function for which /(0) = /(l) = must return 
1 in N/2 of the remaining N — 2 arguments. The number of 
ways in which this arises is (jyTg)- Similar arguments apply 
to the other cases. 



Thus, if x ^ y then 



E c-u 

/ balanced 



N - 2\ / N - 2 
N/2 ) ~ \N/2 - 1 

iV-2\ 2 
N/2 J N/2 - 1 
B 



N - 1 



(7) 



where B is the number of balanced functions and the last 
two lines follow from algebraic manipulations of combi- 
natorials. Eqs. (J6j) and ([7]) imply that, if x ^ y, 



1 



Pf bal s 



N — 1 



PO xy 



(8) 



The cases of all values of x and y are then summarized 
as 



f 



N-l 



Pf bal 



N - 1 



-po 



N E PxpoPx 



(9) 



x=0 



where P x := |x)(x|. 

Minimum error discrimination between the two density 
operators Pf cons t and Pfbai requires a POVM with two 
outcomes and two positive operator elements 7r cons t and 
7Tbai, where 7r C onst + 7Tbai = I- The probability with which 
an incorrect inference is made is 



Po 



= Pconst Tr [f)f balTTconst] + Pbal Tr [/5f const ^bal] (10) 



where p C onst (pbai) is the probability of selecting a func- 
tion from the constant (balanced) class. A standard 
derivation [l2j gives 

Pcrror = jr - 1 1 Pconst Pf const — PbalPf bal||) (H) 

where the trace norm satisfies 



i|| :=Trk/itI] =E^(^) 



(12) 



with {CTi(^4)} being the singular values of A. The trace 
norm is clearly invariant under unitary transformations 
in the sense that, if V is any unitary then \\A\\ = 
||yAV rt ||. Although Eqs. (JTHJ and (JTTJ> are equivalent, 
they have distinct uses in terms of determining condi- 
tions under which the algorithm will succeed. As we 
shall show, Eq. ((TTj) yields the optimal initial state and 
Eq. (fT0|) . the optimal measurement for success. 

Eq. ((TTj) implies that the quantum algorithm for solv- 
ing the Deutsch-Jozsa problem will succeed with cer- 
tainty when ||pconstPf const - PbalPf bal 1 1 = 1- Note that 

llPconstPf const || = Pconst and ||p b alPf bal || = Pbal, giving 
llPconstPf const || + HPbalPf bal || = Pbal + Pconst = 1- Thus 

the quantum algorithm will succeed with certainty if and 
only if 



iPconstPf const ~ PbalPf ball 



I Pconst Pf const || + IjPbalPf ball 



(13) 
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In this context, an important general result [13] is that 
if A and B arc positive semidefinite operators then 
||i - B\\ = ||i|| + ||B|| if and only if A and B have or- 
thogonal support. The support of a positive semidefinite 
operator is the subspace spanned by its eigenstates which 
correspond to non-zero eigenvalues. Thus the quantum 
algorithm will succeed with certainty if and only if the 
operators Pf const and /5f Da i have orthogonal support, or 
cquivalently 



Pi const /3f bal = PfbalPf const — 0. 



(14) 



Defining A := ^2^=0 PxPqPxi which is easily shown to 
be a positive operator, and using Eqs (J3)), (J9)) and (|T4"|) 

gives 



[po, A] = and 
NAp = pi 



(15) 
(16) 



Eq (j!5[l implies that po and A can be simultaneously 
diagonalized. Denote the associated basis of eigenstates 
l>.v{o. j 1 V|, Thus 



P<3 = E 7 '^'^' 

3=1 



(17) 



where L > is the number of non-zero eigenvalues of po 
and < rj ^ 1 satisfy X^t=i r J = Likewise 



■13- 

N 



A = X! A 3^3)(^' 
3=1 



(18) 



where Aj ^ 0. Eq. (|T6|) implies that NXjVj = r| for 
j = 1, . . . , L and this gives 



for j = 1, . . . , L. Additionally, 

N-l L 

A = E^E r ^K^ 

x=0 j=l 
N-l L 

= E^E*^i 2 

x=Q j=l 

where 4>j(x) := (x\<f>j) . Thus Eq. (fT8|) gives 

N-l L 

^ = EE^i^wi 2 i^^)i 2 

x=0 j=l 

and, when combined with Eq. (fl9|) . 

JV-l L 



(19) 

(20) 
(21) 



(22) 



77 = EE*-wi 2 i^)i' 

x=0 j=l 



(23) 



for k ^ L. The fact that L > implies that n 7^ 0. Thus 

N-l N-l L 



jj=riJ2 \M*)\ 4 + E E *(x)|%(x)| 2 (24) 

x=0 x=0 j=2 



where the second term on the right is defined to mean 
when L = 1. The second term on the right of 
Eq. (|24|) is non-negative. The first term contains 
X^ r =o 1 |0i(x)| 4 which is subject to the constraint that 
12^=o \4>i( x )\ 2 = 1- A Lagrange multiplier approach 
shows that Y^x^o^i 00 )^ ^ V-^; provided that < 
|<^i(^)| 2 ^ 1) which is always satisfied. This mini- 
mum is attained when |(^i(x)| 2 = 1/N (or cquivalently 
<j>x{x) = e w * /VN where 9 X is real) for x = 0, . . . N - 1. 
Thus the second term on the right of Eq. is identi- 
cally zero and L = 1. Thus the quantum algorithm will 
succeed with certainty if and only if the input is a pure 
state, i.e. (dropping the subscript and changing notation 
to be consistent with Fig. [T]) 



po = |* )(*o| 



(25) 



The pure state must be such that it produces a minimum 
for the first term on the right of Eq. (|24|) and thus 



I*, 



N-l 
x=0 



(26) 



Eq. (JTOJ) and the positivity of the POVM elements 
and density operators requires that for the algorithm 
to succeed with certainty the measurement operators 
must satisfy Tr [p f bai^onst] = Tr [/5 f const 7r b ai] = 0. But 
TTconst = I — 7Tbai implies implies that this requirement is 
equivalent to Tr [p f const 7r cons t] = Tr [p { bai^bal] = 1- The 
fact that the density operators and the measurement op- 
erators are positive semidefinite and that their eigenval- 
ues fall in the range [0, 1] then implies that 7T C onst must be 
the projector onto the support of /if const- The associated 
POVM elements are 



TTconst = |*o) (*o I and 

TTbal =1- |*o) (*0|- 



(27a) 
(27b) 



Eqs. (|2"6"j) and (|2T|) give the initial states and measure- 
ments for a general quantum algorithm which solves the 
Deutsch-Jozsa problem using a single invocation of an 
oracle unitary of Eq. (JTJ) and no ancillary qubits. The 
standard algorithm, for which |^ ) = -jj^J2x=o\ x )' * s 
one example of this [TH, [l6| . Also the general algorithm 
will successfully identify the function class regardless of 
the probabilities with which the various admissible func- 
tions are selected. This can be verified by applying the 
algorithm unitary transformations for the various admis- 
sible functions to the initial state of Eq. (f2"6")) , computing 
the final state and using this to determine the POVM 
operators of Eq. ([27]) to determine probabilities of the 
two possible outcomes. 
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IV. APPLICATION TO THE DEUTSCH-JOZSA 
ALGORITHM ON MIXED INITIAL STATES 

In some proposed implementations of quantum com- 
puting, such as solution state nuclear magnetic resonance 
(NMR) [l8r:27|, the initial state of the quantum system 
is mixed and therefore the problem cannot be solved with 
certainty by using the scheme of Fig. [TJ The notion of 
minimum error discrimination between quantum opera- 
tions can be applied to bound the success probability of 
any algorithm on such mixed input states. For a general 
mixed state, Eqs. © and ((TTJ) imply 



1 - 



Pconst 



Pbal \ * PbalN ~ 

Po - — — t A 



N - 1 



N - 1 



(28) 

We consider the special case where p C onst = Pbai = 1/2 
and thus 



N 



2(N- 1) 



Po- A 



(29) 



Attempting to apply the triangle inequality to the trace 



A 



1 results in 



norm and using the facts that \\po\\ 

an inequality which is always satisfied. 

However, in some circumstances a similar approach can 
yield meaningful bounds. Specifically we consider situa- 
tions in which an ensemble of quantum systems is first 
allowed to reach thermal equilibrium followed by appli- 
cation of a preparatory, oracle-independent unitary. The 
resulting state constitutes the initial state for the system. 
We shall assume that the system consists of n qubits and 
has a Hamiltonian 



H 



E 



9 z 



(30) 



where w, is the precession frequency of the i th qubit and 
Jij is the coupling between the i th a nd 7 th qubits (this 
is typical for solution state NMR dlGl]). Tlie tnerrnal 
equilibrium density operator is pth = e~^ H /Z where j3 = 



1/fcTandZ = Tr 



-0H 



We shall assume that J, 



'ij < u k 

and we shall ignore coupling terms in the density operator 
which is constructed from 



-ftH 



(31) 



where on := hui/2kT. Our final assumption is that aj -C 
1 (in typical solution state NMR scenarios, u>i ~ 500 MHz 
and T ~ 300 K so that a, ~ 10" 5 ). Thus 



Pth 



1 - 
N 1 



Pdcv 



(32) 



where the deviation density operator is pdcv '■= 
J2i=i ati&z /N. Note that Tr/Sd cv = 0. Given the thermal 
equilibrium state as an initial state of the system, there 
is a range of possible quantum operations that can be 



applied prior to the first invocation of the oracle. These 
could contain non- unitary operations; examples include 
pse udop ure state preparation [2 Ell or algorithmic cool- 
ing |28| . Our aim is to focus on scenarios in which such 
non-unitary operations are avoided and we only consider 
scenarios where operations applied prior to the first ora- 
cle invocation are unitary without involving any ancillary 
systems. 

Here the most general initial state of the n qubit sys- 
tem is represented by p = Vp t hV^ where V is the 
preparatory unitary operation. Then 



1 



Po 



where Pdev 



N 

Vp dcv V^ and 

N-: 

f - 

N 



Pdc 



A 



— E PxPdcvPx- 



(33) 



(34) 



Define A dev 



Pc 



x=0 

I'Pd.J',- Thus 
N 



1 



2(N- 1) 
The triangle inequality yields 



Pdc 



A 



dev 



1 



f 



N 



2(N- 1) 
The diagonal nature of A dov gives 









(llpdcvll + 


A' 


)] 



(35) 



(36) 



N-l 



A' 



(lev 



Ek^c 



Pdcv I x ) 



(37) 



Denote the eigenstates and eigenvalues of p' dcv by 
{|Xt} I i = 1, . . . , JV} and fa \ i = 1, . . . , N} respectively. 
Then 



A' 



dev 



N-l N 

\ci\\(x\Xi)Y 

x=0 i=0 
N N-l 



EE 

x=0 i=0 
N A 

i=0 x=0 
N 

Ei 



I Pdc 



Thus 



Pc 



1 - 



N 



(N-l) 



I Pdc 



(38) 



(39) 



The unitary invariance of the trace norm 
implies that ||p dcv | = ||pdcv|| • The singu- 
lar values of pdcv are {\a± + a 2 + . . . a n \/N, 
\—ai + a 2 + ... a n \/N, |ai - a 2 + . . . a n \/N, . . .} . 
Without loss of generality assume that ol\ ^ a 2 ^ 
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o;3, ■ • ■ ^ 0. Then each singular value is bounded from 
below by na\/N and 



Thus 



IPdevli < noti. 



1 / Nnax 

Perror> £ I 1 ~ 



(40) 



(41) 



Defining e := Nna±/(N — 1) gives p e rror ^ (1 — e)/2. 
This gives the failure probability for an application of 
the algorithm to a single ensemble member and it indi- 
cates a non-deterministic output. The algorithm may be 
run repeatedly using the same unitary on a large number 
of independent quantum systems, all initially described 
by the same density operator and this will ultimately 
increase the chances of successfully identifying the ora- 
cle class. However, this must be compared to a classical 
probabilistic algorithm. Such an analysis has been done 
in the case of a pseudopure-state input [23 ]. Here the 
quantum algorithm is reconfigured, using a single ancil- 
lary qubit to which the function class is written when- 
ever a correct pure state input is used. For mixed input 
states, it is possible that an erroneous function class may 
be written to this ancillary qubit. Over an entire en- 
semble, the function class is inferred by effectively taking 
a "majority- vote" of computational basis measurement 
outcomes on the ancillary qubits for individual ensemble 
members. The probability of successful inference only 
depends on the ensemble size and the probabilities with 
which each of the two possible measurement outcomes 
occur; the initial state merely determines these probabil- 
ities in terms of e. This is to be compared with a classical 
probabilistic algorithm. The results are that, for large en- 
semble size (i.e. in the limit as this becomes infinite) the 
quantum algorithm succeeds with larger probability than 
the classical probabilistic algorithm if e > \/3/4. For the 
quantum algorithm to succeed with greater probability 
than the classical algorithm, this implies 



ai > 



N-l 1 



N 



or, alternatively 



n > 



N 



N 



1 1 

ai 



In typical solution-state NMR situations a\ 
thus 



n > 



N-l 



N 



10 5 



(42) 

(43) 
10" 5 and 

(44) 



Thus with current solution state NMR technology and 
an implementation of the algorithm which starts with 
the thermal equilibrium state, uses no ancillary qubits 
and only uses the oracle unitary of Eq. ([1]) plus any 



other unitaries, the minimum number of qubits required 
for the quantum algorithm to succeed with greater prob- 
ability than any classical algorithm is approximately 
10 5 . Current implementations of the Deutsch-Jozsa 
algorithm have not even exceeded n = 10. 

V. CONCLUSION 

In conclusion we have demonstrated the usefulness of 
the notion that quantum algorithms are tools for discrim- 
inating between oracle unitary transformations. For the 
Deutsch-Jozsa problem, applying techniques associated 
with unitary discrimination results in the complete set 
of quantum algorithms which can solve the problem with 
certainty; these have a structure similar to that of the 
standard version of the algorithm. We have also explored 
the issue of solving the Deutsch-Jozsa problem when the 
initial state of the quantum system is selected from a re- 
stricted set. For the case of a thermal equilibrium state 
followed by an arbitrary oracle-independent unitary, this 
yields a lower bound on the number of function argu- 
ments such that the quantum algorithm succeeds with 
greater probability than any classical algorithm. This is 
several orders of magnitude larger than that which has 
been implemented experimentally to date. 

We should point out that our analysis assumed cer- 
tain constraints. First wc assumed the oracle unitary 
of Eq. ([T]), which operates on n qubits. However, some 
versions of the algorithm require n argument qubits plus 
one ancilla qubit and use an oracle unitary defined by 
Uf\x)\y) := \x)\y © f{x)) (the rightmost system in this 
notation is the single ancilla qubit, the leftmost the n 
argument qubits) [7|. The unitary considered in our 
analysis can be reached from this by fixing the ancilla 
qubit in a special state, converting the problem to one of 
phase estimation [ll| . Whether this yields a larger set of 
quantum algorithms which solves the problem with cer- 
tainty or reduces the lower bound in the thermal equilib- 
rium state scenario is open to investigation. Second, it is 
known that using ancilla qubits entangled with "system" 
qubits enhances the possibility of successful discrimina- 
tion between unitaries [llf and it is conceivable that this 
could also yield a larger set of quantum algorithms and 
improved lower bounds in the thermal equilibrium case. 
Finally, in the thermal equilibrium scenario, relaxation 
to the thermal equilibrium state could be followed by a 
non-unitary quantum operation (30j and we have not as- 
sessed bounds in any such cases. 
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